Andreev reflection and cyclotron motion at superconductor - normal-metal interfaces 
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We investigate Andreev reflection at the interface between a superconductor and a two-dimensional elec- 
tron system (2DES) in an external magnetic field such that cyclotron motion is important in the latter. A finite 
Zeeman splitting in the 2DES and the presence of diamagnetic screening currents in the superconductor are 
incorporated into a microscopic theory of Andreev edge states, which is based on the Bogoliubov-de Gennes 
formalism. The Andreev-reflection contribution to the interface conductance is calculated. The effect of Zee- 
man splitting is most visible as a double-step feature in the conductance through clean interfaces. Due to a 
screening current, conductance steps are shifted to larger filling factors and the formation of Andreev edge 
states is suppressed below a critical filling factor. 
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I. INTRODUCTION 

The trend toward radical miniaturization of electronic de- 
vices has spurred efforts aimed at understanding the effect of 
quantum coherence on transport. Current research is inves- 
tigating a variety of such mesoscopic systemsii 2 ^ Recently, 
the subject of mesoscopic superconductivity^^ has attracted 
considerable interest, partly due to the possibility to perform 
experiments in hybrid structures consisting of superconduc- 
tors and semiconductors where phase-coherence lengths ex- 
ceed device dimensions. 7 Many effects observed in these sys- 
tems ulimately stem from Andreev reflection?- 9 Speaking in 
simple terms, Andreev reflection occurs at the interface be- 
tween a superconductor (S) and a normal metal (N) when an 
electron with energy within the superconducting gap is inci- 
dent on the interface from the N side. While no states are 
available in S for the incoming electron, a finite amplitude 
for pairing with some suitable electron from N is induced by 
the proximity of S, which allows both electrons to enter S as 
a Cooper pair. If that happens, N is left with a hole excita- 
tion that has the sign of its group velocity, charge, and mass 
opposite to that of the electron. As incoming electron and 
Andreev-reflected hole have a definite phase relation, bound 
states are formed in multi-interface geometries such as S-N- 
Si£ or S-N-tii systems, giving rise to discrete energy levels 
within the superconducting gap. 

An external magnetic field B = V x A affects transport 
properties of mesoscopic systems at least in two ways. In sys- 
tems with finite Zeeman splitting (given by gftsB, g being 
the gyromagnetic factor and [is the Bohr magneton), spin de- 
generacy is lifted, which turns out to suppress Andreev reflec- 
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In addition, the vector potential A explicitly enters 



the single-electron Hamiltonian, 

#o = ^- (p - elf + | 3 ■ B 
2m V / 2 



(1) 



resulting in a host of peculiar quantum-mechanical effects. 
One of the most well-known is perhaps the appearance of 
phase factors which lead to magnetoconductance oscillations 
in a variety of mesoscopic systems 3 and affect orbital wave 



functions even in regions where no magnetic field is present. 15 
If the magnetic field is treated within a quasiclassical approx- 
imation which assumes the relevant physical length scale of 
the system to be much smaller than the cyclotron radius of 
electrons, the Hamiltonian of Eq. Q reduces to 



H^ c> = e F + v F p ■ (p-eA 



fi B a ■ B 



(2) 



and modifications of orbital wave functions due to the mag- 
netic field arise entirely from such phase factors. Here, e F = 
h 2 kp I (2m), up = hkp/m, and k F are the Fermi energy, Fermi 
velocity, and Fermi wavevector, respectively. Most of the pre- 
vious studies of Andreev reflection in a magnetic field were 
carried out in the quasiclassical HmitiAiii&Aii&iSiSSi^ii^SiS 3 . 

The quasiclassical approximation ceases to be valid, how- 
ever, when the cyclotron radius r c = l B k F is smaller than the 
characteristic length scales of the system which are set, e.g., 

-* 1/9 

by disorder, temperature, or size (here, l B = (h/\eB\) ' de- 
notes the magnetic length which is the fundamental quantum- 
mechanical length scale introduced by the magnetic field.) In 
this limit, another well-known quantum effect due to the mag- 
netic field becomes important, which is Landau quantization 24 
of electron motion in planes perpendicular to the field direc- 
tion. In transport, Landau quantization manifests itself, e.g., 
by Shubnikov - de Haas oscillations and by the quantum Hall 
effect. 25 ' 26 Recent advances in fabrication technology make 
it possible to study Andreev reflection in superconductor- 
semiconductor hybrid systems at magnetic fields where cy- 
clotron motion of quasiparticles can be important in the semi- 
conducting regionpi^LSi^ii^i motivating a number of related 
theoretical studies, 32 . 33 ! 34 . 35 . 36 ! 37 . 38 . 39 ! 40 

Here we present results on the interplay between Andreev 
reflection and cyclotron motion at S-N interfaces. We focus 
in particular on the effect of a finite Zeeman splitting in the 
normal region and on the impact of the diamagnetic screening 
current in the superconductor.— To be specific, we choose a 
geometry where B = B z and the interface is the yz plane. 
Considering particle motion in the xy plane from a classical 
point of view, the interplay between cyclotron motion and An- 
dreev reflection at an ideal interface results in electrons and 
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holes alternating in skipping orbits along the interface 42 ^ in y 
direction, as shown in Fig. [2 In quantum mechanics, however, 
the guiding-center coordinates for cyclotron orbits are canon- 
ically conjugate operators which cannot be diagonalized si- 
multaneously 44 while each of them separately commutes with 
Hq. In the present problem, where translational invariance is 
broken due to the interface, it is useful to choose the represen- 
tation where the guiding-center coordinate X in x direction is 
a good quantum number, i.e., eigenstates of Ho are labeled by 
X. Details of our quantum-mechanical description of the S-N 
interface in a magnetic field are given in the following section. 
It is based on finding solutions of the Bogoliubov - de Gennes 
(BdG) equation^ 



H {) + U- e F 



IAI e 2 ^ 



lAle" 2 ^ -m-U 



= E 



(3) 

for the hybrid system by matching at the interface appropriate 
solutions for the S and N regions. 46 Scattering at the interface 
is modeled by an external potential U = U{x) = Uq S(x). We 
allow for different effective masses and Fermi energies in the 
S and N regions. For our study, we adopt the model where 
the relevant quasiparticle excitations in the S region are those 
of a conventional BCS s-wave superconductor with a uniform 
superflow in the y direction parallel to the interface due to the 
presence of diamagnetic screening currents. This approach is 
strictly valid only when S is a type-II superconductor in the 
Meissner or diluted-vortex phases. To be specific, we con- 
sider exactly this case in Sec. |II] but discuss generalizations 
to the mixed phase in the Discussion section. The motion in 
z direction, i.e., parallel to the magnetic field, can be trivially 
separated and gives only rise to an additive renormalization of 
the Fermi energy. We shall neglect this in the following, mak- 
ing our study especially applicable for S-N junctions where 
the N region is a two-dimensional electron system. Eigenen- 
ergies for the motion in the xy plane in a uniform magnetic 
field are labelled by the guiding-center quantum number X 
and correspond to the familiar Landau levels of independent 
electrons and holes when \X\ ^> r c . Similar to the case of 
a hard wall, 47,48 Landau levels are bent for X close to the S- 
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FIG. 1: Semiclassical picture of Andreev edge states. Electrons and 
holes in skipping orbits along a clean S-N interface are successively 
Andreev-reflected into each other. In the S region, quasiparticle ex- 
citations exist only on the (extremely short) length scale of the su- 
perconducting coherence length. The presence of a uniform super- 
current parallel to the interface results in a shift of electron and hole 
orbits in the normal region. Andreev reflection couples electron and 
hole orbits with guiding-center coordinate ±X — A, respectively. A 
denotes the magnetic penetration depth. 



N interface, and the corresponding eigenspinors of the BdG 
equation are of mixed electron-hole type. An example of the 
spectrum of such solutions with energies within the supercon- 
ducting gap was given in Ref.l33l 

The classical picture of skipping orbits along the inter- 
face, and the analogy with the familiar quantum-Hall edge 
states, 47 - 48 - 49,50 suggests that currents will flow along the in- 
terface in the plane perpendicular to B (i.e., in y direction). 
In superconductivity, quasiparticle current j P = j e + j\ and 
charge current jq = e(j e — jh) are distinguished where, for 
a particular solution of the BdG equation, electron and hole 
currents are given byi 5 - 

J e = — 3?e ju* (p - eAj uj , (4a) 
% = —^fte^v (p-eAj w*} . (4b) 

It turns ou*2£ that j e and jh are always parallel to the inter- 
face for states with definite guiding-center coordinate X in 
the direction perpendicular to the interface. This is no surprise 
because, owing to the peculiar quantum dynamics in a mag- 
netic field, the guiding center for electrons and holes cannot 
be localized in y direction in a state where it is fixed in the x 
direction. Generalizing the Hellmann-Feynman theorem 5 ^ to 
the matrix BdG Hamiltonian of Eq. Q, a general expression 
can be derived 33 for the contribution of a particular eigenstate 
(unX , v n x) with energy E n x, labelled by Landau-level index 
n and guiding-center quantum number X, to the total quasi- 
particle and charge currents in y direction: 



r(P,Q) 
l nX 



(jP,C 



(5) 



The total quasiparticle current turns out not to depend explic- 
itly on the electron or hole wavefunctions. Rather, it is given 
entirely in terms of the energy spectrum, 



r(P) _ - _ 
1 nX ~ • - 



h 1 



8E r 



X 



dX 



(6) 



which aquires a nonzero dispersion only close to the S-N in- 
terface. Hence, confinement due to the interface with a super- 
conductor gives rise to the formation of a chiral edge channel 
for quasiparticle current in exactly the same way as confine- 
ment due to a hard wall induces edge channels for charge cur- 
rent^ 4 ^ We distinguish the new kind of edge channel real- 
ized at the S-N interface from its counterpart near a hard wall 
by calling it an Andreev edge channel and call the solutions 
(u n x, v n x) of the BdG equation with nonzero Irt Andreev 
edge states. [In Eq. (|6j, L y denotes the system size in y di- 
rection.] As expected intuitively, the total charge current I^ x 
flowing in an Andreev edge state does depend on the details 
of the wave functions; it can be written as I^ x = ffiL + I^ x . 
The expression for I®r x suggests a straightforward interpreta- 

33 



tion. It is given by 



t-(H) 

1 nX 



dE r 



X 



h L„ dX 



(1-25, 



nX 



(7) 
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where B n x — J \v n x\ 2 is the hole probability in the An- 
dreev edge state. Equation is the expression for current 
in an edge channel that is reduced due to the presence of An- 
dreev reflection. For vanishing B n x, i e, when only normal 
reflection occurs at the S-N interface, Andreev edge chan- 
nels actually become the familiar quantum-Hall edge chan- 
ne j Si 47 j 48 i 49 i 50 j nS p ec tion shows that the term represents 
the quasiparticle-conversion current in the superconductor 
that supports the Andreev edge state. Note that the above 
discussion, as well as the analytic expressions, for currents 
are generally valid, i.e., in particular also with finite Zeeman 
splitting and diamagnetic screening currents present. It is then 
possible^ to express the Andreev-reflection contribution to 
the S-N interface conductance in terms of the hole probabil- 
ties of zero-energy Andreev edge states as 



G 



ar — — rT^ B n , 



(8) 



where the label n comprises the set of quantum num- 
bers (Landau-level index and spin-projection eigenvalue) 
for zero-energy eigenstates, and the indicates that only 
electron-like states (i.e., those with B n < 1/2) are to be 
summed over. As the experimentally relevant quantity, we 
present results for Gar in the following. 

We continue with presenting our model of an S-N junc- 
tion subject to a quantizing magnetic field in Sec. |ll| (Read- 
ers not interested in the mathematical formalism can skip this 
section.) Numerical results obtained within that model are re- 
ported in Sec. [HI] where the effects of Zeeman splitting and 
diamagnetic supercurrents are elucidated in detail. We dis- 
cuss implications of these results with particular emphasis on 
experimentally realizable superconductor-semiconductor hy- 
brid structures. Conclusions are given in the final Sec. lIVI 



II. MODEL OF AN S-N HYBRID SYSTEM IN A 
QUANTIZING MAGNETIC FIELD 

Our model of the S-N interface is similar, in spirit, to that 
used in the approach by Blonder, Tinkham, and Klapwijk 46 
(BTK) and its recent generalizations to superconductor- 
semiconductor interfaces. 52 We consider a planar hybrid sys- 
tem (in the xy plane) consisting of a semi-infinite two- 
dimensional electron system (2DES, located in the half-plane 
x > 0) and a superconductor (occupying the half -plane x < 
0). The effective mass of electrons, the Fermi energy, and the 
modulus of the superconducting pair potential are assumed to 
be piecewise constant as functions of the coordinate perpen- 
dicular to the interface: 

m{x) = m s 0(— x) + m n Q(x) , (9a) 
e F (x) = e F>s 6(-a;)+e F)n B(x) , (9b) 
|A|(ar) = A 9(-z) . (9c) 

The magnetic field, applied perpendicular to the 2DES, is as- 
sumed to be screened from the S region, 



B(x) = {B Q(-x) + B B(x)} z 



(10a) 



where A denotes the magnetic penetration depth. Transla- 
tional invariance in y direction suggests a gauge for the vector 
potential where A = A(x) y, and we choose 



A(x) = A B (e* - l) Q(-x) + x B 6(x) 



(10b) 



In the presence of a magnetic field, the phase 2ip of the super- 
conducting pair potential aquires a non-trivial dependence on 
spatial coordinates, and care has to be taken to determine it 
correctly. 1619 Combining Maxwell and London equations, 



V x B 

Js 



MO Js 

h 



(lla) 
(lib) 



\2 s ' 

where J s is the screening supercurrent, and 

k s = Vip-eA , (12) 
we obtain from Eqs. (llOi for the case under consideration 

v(v) = ~w y s s n ( eB ) + vo , (i3) 

L B 

with sgn(x) denoting the sign function, and (po a constant. 

We calculate electronic and transport properties from solu- 
tions of the BdG equation 0. Upon making the separation 
Ansatze 



u(x,y) 



v(x,y) = 



1 



fx{x) e 



» My)+?/ ttt sgn(eB) 



gx(x) e 



» -¥>(»)+!/ Tf-sgn(e-B) 



(14a) 



,(14b) 



and using the fact that, in the presence of Zeeman splitting, 
Eq. l|3} exhibits a block-diagonal form 53 in the spin quantum 
number a, we find a onedimensional BdG equation, 



H+.a - 

IAI 



U |A| 



fx,o 
9X-CT 



E 



fx,o 

9x- 



(15) 

The single-particle Hamiltonians H±, a are different in the S 
and N regions, H±. a = ®(-x) Hf a + Q(x) H^ a , where 



H 



(n) 

±,<T 




(16) 

Here lo c = h/(m n l 2 B ) denotes the cyclotron frequency in the 
N region, rjz — g/^sB / (huj c ) measures Zeeman splitting, and 
v = 2eF >n / '(Hu>c) is the filling factor. In the superconductor, the 
single-particle Hamiltonian depends on the superflow wave 
vector fc s = fc s y, 



A 



sgn(ei?) 6(— x) 



B 



rr(s) 
U ±,a 



jr. t h 2 [X±l 2 B k sS gn(eB)Y 



2m 



2m s l% 



(17a) 



e F ,s<17b) 
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Here we have neglected Zeeman splitting in the S region, 
which is a good approximation for typical field strengths ap- 
plied experimentally in S-2DES hybrid structures. For our 
purposes, we need to find solutions of Eq. (I15> in the S region 
on length scales x ~ £o where £n is the superconducting co- 
herence length. In a type-II material which has £o *C A, we 
can approximate expjx/A} w 1 in Eq. < !17al and use 



rr(s) 

a ±,o 



2m s 



2m s lt 



(X T A) - e F>s 



(17c) 



Our approximation amounts to neglecting the cyclotron mo- 
tion of evanescent quasiparticles in the S region, which trans- 
lates into the condition that the superconducting gap parame- 
ter Ao is much larger than the cyclotron energy of electrons 
in the superconducting material. This is typically satisfied for 
realistic situations. 

Solutions of Eq. (1151 for the S-N system are found by 
matching Ansatze for eigenfunctions in the S and N regions 
at the interface x = 0. In the S region, where a uniform su- 
perflow with velocity hX/{m s l 2 B ) is present, the appropriate 
Ansatz is a superposition of Doppler-shifted quasiparticle ex- 
citations, 



x<0 



with functions ij)±(x) = exp{^fixq±}, and parameters 



(18) 



For \E + 5\\ < Ao, the wave function displayed in Eq. Jl 8I > 
corresponds to evanescent quasiparticle excitations in S. The 
Ansatz for the wave function in the N region is a superposition 
of purely electron and hole-like BdG spinors, 



(.^L-(s)^» + (?)^«-»- 

(20) 

where (± = x + X^f X, and the functions x±.&(() are eigen- 
functions of the harmonic-oscillator Schrodinger equation, 



1% d 2 X±, 



' B 

2 



r a2 



d( 2 



C _ v - cryz E 

2 ~ o T 



2l 2 



fiUJc 



X ±.a=0 , (21) 



normalized to unity in the half-space occupied by the N re- 
gion. They can be expressed in terms of parabolic cylinder 
function}: 54 that are well-behaved as x — > oo, 



r i' I v-ariz E R C± 
X±AC±)=F±,.U[ — T —V2- 



(22) 



q± = 
7± = 

5\ = hu>, 



2m s 



e F , s ±^/(E + S x ) 2 -A 2 

An 



E + S x T x/iE + S^-A 2 
m n XX 
m s 1% 



/J 



(19b) 



(19c) 



where F± t<T are normalization constants. 

Imposing the continuity of BdG quasiparticle wave func- 
tion and conservation of the probability current at the interface 
yields the secular equation 



GH(c' 2 + d' 2 ) + G'H' = c' ■ E + Sx ^=(G'H -GH')+d'(G'H + GH') , (23) 

\/A 2 ~(E + S x ) 2 



I 

whose parameters are: d — ^ Sfte {q+}, d' = — 3m {</+}+ normal region: 



U with U=^U ,G = x+,<r(A-X), H = X -,-a(A+X), 
G' = d x x+,o{x + X — X) \ x -o, and H' = d x X-,-a{x + A + 

X) \ x —q. Equation fl23i is an implicit relation between the G' + G I E+Sx c' — d! 



guiding-center coordinate and the excitation energy E, yield- — = y-\/ A o (E+ s \) ^ 

ing the dispersion relation _E vs. X. a c '-^ ^a§-(b+5 a ) 2 



Since it is important for finding the Andreev-reflection con- 
tribution to the conductance, we also give the general expres- 
sion for the amplitude ratio aj b of the wave function in the 



Further analytical progress can be made using the asymp- 
totic expansion 54 for parabolic cylinder functions U(e, x) « 

cos (e/2 + l/4)7ry / |e|a; , being T the Gamma 

function. In our case this expansion is valid when e Fn ± (E — 
ag^B/2) ^> m a uj 2 (\X\ + A) 2 . With this, the approximated 
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secular equation reads 



cos(<p +)0 -)+fi(^_ )0 .) 



25 (E + 5 X ) sm(<p+. a ) 



8.9 



(25) 



i.e., as- 



Here we employed the Andreev approximation 
sumed E, An, \g^B\ -C £F, n , £f,s> and used the abbreviations 



The parameter s = [e^ s m n /(eF, n 'm s )] 1 ^ 2 measures the 
Fermi- velocity mismatch for the junction, and w = 
[2m n UQ /(f^e^n)] 1 / 2 quantifies interface scattering. Within 
the same approximation, we can write the amplitude ratio as 



<P-,<T 

Cl(a) 



E _ arjz 
Hcu c 2 
2 



E 



2jD 
2 



X 
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' E 








2 



r -2V^ r 



+ w 2 — 1] sin(a) + 2w cos(a) 



(26b) 
(26c) 



1 - sin(y +)(T + 



- =J\f- 



E+S x 



[sm((p +t<7 ) - cos(^_ )CT )] 



(27) 



where 



^ = r[l/4(l + v - a m ) + g/(2guQ] 2E/(toc) 



r[i/4(i 



In the following, we discuss the qualitative effects of fi- 
nite Zeeman splitting and screening supercurrents separately, 
based on our approximate analytic approach described above. 
While this allows us to understand certain basic features ex- 
hibited in the exact numerical results to be presented later, we 
caution the reader that the quantitative estimates given in the 
remainder of this section can be expected to apply only in lim- 
iting cases that are consistent with our approximations. 



A. Finite Zeeman splitting, no screening current 

Here we discuss briefly the effect of the Zeeman splitting 
on Andreev edge states in the absence of screening currents. 
We consider the case of small excitation energies, i.e. \E — 
ag^s,B\ <C hy / e^ n /(2m n X 2 ). In this situation, we find 



E a = A 



(2n + 1)tt ± arccos(Slo) — 2X^/2m n e^ n /h — 



q + 7rAo/ (Hlu c ) 



where On = f2(?rz//2) and q = 2s/(s 2 



w 



1) 



(28) 
Equa- 



tion A28I is very similar to the results obtained in the absence 
of Zeeman interaction in Ref. |33[ In particular, it shows that, 
at low excitation energy (those that are relevant for linear 
transport), the Zeeman splitting can be absorbed in a spin- 
dependent shift of the guiding-center coordinate X. We there- 
fore expect the Zeeman splitting to result in the usual halving 
of conductance steps as a function of filling factor. As we 



will see below, however, this feature manifests itself only for 
very transparent interfaces, while Zeeman splitting turns out 
to have a small effect on the conductance in realistic situa- 
tions. 



B. Effect of the screening current 

It is immediately apparent from the general form of the sec- 
ular equation that there are two important consequences of fi- 
nite screening currents. First, a shift of quasiparticle energies 
in the superconductor by 5\ occurs. [See Eq. ( I19ci .l As An- 
dreev edge states have only evanescent BdG spinor wave am- 
plitudes in the superconductor, their formation is now possible 
only when the condition Ao > \E + 5\ \ is satisfied. Second, 
the spatial position corresponding to zero guiding-center co- 
ordinate of coupled electron and hole wave functions in the 
2DES is shifted to a; = —A into the superconducting region 
(see Fig. 0. Since the matching of BdG spinors still must be 
performed at the S-N interface (x = 0), Andreev-reflection 
coupling in the normal regions becomes less direct. 

For a more quantitative assessment of the effect of screen- 
ing currents on the Andreev-reflection conductance, we start 
by discussing the ideal interface, i.e. w — and s = 1. In this 
case the secular equation reduces to 



cot(y_| 



E + 5 



A 



(29) 



For zero-energy bound states (E = 0) and in the absence of 
Zeeman splitting, Eq. ( I29i can be transformed into the secular 
equation for the system without screening current but at finite 
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energy S\ and with a renormalized filling factor 



7T m n Afc F>n 

v I 1- 

2 m= v 



(30) 



This correspondence implies that sharp features in the linear 
conductance as a function of filling factor, which have their 
origin in the sudden appearance of new bound-state levels, are 
shifted to larger filling factors due to the presence of a finite 
screening current. In addition, the formation of zero-energy 
Andreev edge states will be suppressed below a critical filling 
factor 



m n 6 F ,n 

m s A 



Afc F 



(31) 



because then their quasiparticle amplitudes in the supercon- 
ducting region will not be evanescent. Hence, the Andreev- 
reflection contribution to the interface conductance vanishes 
for v < v„. 



In the more typical case of a nonideal S-N interface, i.e., 
one with a nonzero interface resistance, the screening current 
also renormalizes the effective interface parameters s and w. 
Considering again zero-energy edge states in the absence of 
Zeeman splitting, we find new parameters 



2wcos(2fc Fn A) - (s 2 + w 2 - 1) sin(2/c Fn A) 
s 2 + w 2 + 1 - (s 2 + w 2 - 1) cos(2fc Fn A) - 2wsin(2fc Fn A) ' 

2s 

s 2 + w 2 + 1 - (s 2 + w 2 - 1) cos(2fc Fn A) - 2wsin(2fc Fn A) ' 



(32a) 
(32b) 



III. NUMERICAL RESULTS AND DISCUSSION 

We have determined the BdG spinor wave functions for 
zero-energy Andreev edge states numerically and calculated 
the Andreev-reflection contribution to the S-2DES interface 
conductance according to Eq. Here we present results 
which elucidate the effects of finite Zeeman splitting and su- 
perconducting screening currents. 

For the sake of clarity, we start by considering finite Zee- 
man splitting in the absence of screening currents and for an 
ideal interface. As can be seen from the dot-dashed curve 
in Fig. |2 a double-step structure emerges analogous to the 
conductance of spin-resolved quantum-Hall edge channels. 50 
This feature is quickly suppressed, however, as the size of the 
superconducting gap increases, as well as the interface barrier 
becomes more opaque(w > 0). Using parameters that sim- 
ulate a realistic NbN-InGaAs hybrid structure, Fig. [5] shows 
that the effect of Zeeman splitting can be expected to be rather 
marginal in typical samples. It results only in a suppression of 
the peak conductance, while the oscillatory behavior as func- 
tion of filling factor remains unaffected. 

Turning to the investigation of the diamagnetic screening 
current, we show numerical data for an ideal interface in 
Fig. |4] As expected from our analytical results presented in 
Sec. IIIBI we find that step features get shifted to higher fill- 
ing factors as the magnetic penetration depth increases. The 
magnitude of the conductance is only marginally affected. 
Interestingly, an improvement of step quality is clearly ap- 
parent for increasing penetration depth. As in the case with- 
out screening current, step-like features disappear and are re- 
placed by an oscillatory conductance as soon as the interface 
is nonideal. This can be seen in Fig. [5] where the effect of 
interface scattering in the presence of screening currents is 



illustrated. While some of the sharp features of the step struc- 
ture remain, the Andreev-reflection contribution to the con- 
ductance is much suppressed by an interface barrier. It does 
reach a value comparable to the ideal case (w = 0, s = 1) 
only at specific values of filling factor. 

Finally, the behavior of the conductance in a realistic S- 
2DES hybrid structure with screening supercurrents present is 
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FIG. 2: Effect of Zeeman splitting on the Andreev-reflection con- 
tribution to the S-N interface conductance Gar- A double-step 
structure is clearly visible for the dot-dashed curve which has been 
calculated for an ideal interface (s = 1, w = 0) and g — —20, 
Ao = 0.3 meV, ef.s = £F,n = 10 meV, and m s = m„ = 0.035m c . 
This feature turns out to be obscured at larger values of the supercon- 
ducting gap (Ao = 3 meV for all solid curves), in particular when 
an interface barrier is present (nonzero values of w as indicated for 
each curve). 
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shown in Fig. [6] Vanishing of the conductance below a criti- 
cal filling factor arises from the disappearance of zero-energy 
bound states due to the Doppler shift of quasiparticles in the 
superconductor which reduces the excitation gap to zero4£ As 
the quasiparticles become propagating in the S-region, An- 
dreev reflection at the interface is suppressed. This has impor- 
tant consequences for experimental investigation of Andreev 
edge states. The results shown in Fig.[6]indicate that, for a S 
electrode made from a NbN film for which A ~ 400 nm has 
been reported in the literature, 55 effects due to Andreev reflec- 
tion at the interface may be visible only at rather large values 
of the filling factor. 

IV. CONCLUSIONS 

We have investigated the interplay between Andreev re- 
flection and cyclotron motion at S-N interfaces. The effects 
of Zeeman splitting in the normal region and diamagnetic 
screening currents in the superconductor have been taken into 
account. In the ideal case, Zeeman splitting results in a dou- 
bling of step features in the Andreev-reflection contribution 
to the interface conductance, and screening currents lead to a 
shift of such features to larger filling factors. However, in real- 
istic S-2DES hybrid structures, the Zeeman splitting turns out 
to merely suppress peak conductance. Screening currents turn 
out to be important because the Doppler shift of quasiparticles 
in the S-region will suppress their energy gap to zero below 
a critical filling factor, making Andreev reflection disappear 
altogether. Above this critical value, the typical oscillatory 
structure of the Andreev reflection contribution to the inter- 
face conductance is displayed. 
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FIG. 3: Effect of Zeeman splitting for a realistic S-2DEG junction 
realized, e.g., by a typical NbN/InGaAs hybrid system: ep -s = 7 eV, 
£F,n = 10 meV, m s = m c , m n = 0.035m e , Ao = 3 meV, w = 0, 
and for various values of the gyromagnetic factor g as indicated. 
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FIG. 4: Diamagnetic screening currents lead to a shift of conductance 
steps. Data shown are calculated for an ideal interface (s = 1, w = 
0), g — 0, Ao = 0.3 meV, 6f, s = £F,n = 10 meV, and m s — m n = 
0.035m e for various values of the magnetic penetration depth A. 
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FIG. 5: Andreev-reflection transport at high magnetic fields through 
a nonideal S-N interface with diamagnetic screening current present. 
Results shown are for A = 5 nm and various values of w as indicated, 
with all other parameters as given in Fig.|4] 
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